A PROOF OF PARISI'S CONJECTURE ON THE RANDOM 
ASSIGNMENT PROBLEM 



SVANTE LINUSSON AND JOHAN WASTLUND 

Abstract. An assignment problem is the optimization problem of find- 
ing, in an m by n matrix of nonnegative real numbers, k entries, no two 
in the same row or column, such that their sum is minimal. Such an 
optimization problem is called a random assignment problem if the ma- 
trix entries are random variables. We give a formula for the expected 
value of the optimal fc-assignment in a matrix where some of the en- 
tries are zero, and all other entries are independent exponentially dis- 
tributed random variables with mean 1. Thereby we prove the formula 
l-|-l/4-f 1/9 + - ■ ■ + l/k'^ conjectured by G. Parisi for the case k — m — n, 
and the generalized conjecture of D. Coppersmith and G. B. Sorkin for 
arbitrary k, m and n. 



1. Introduction 

The problem of minimizing the sum of k elements in a matrix of non- 
negative real numbers under the condition that no two of them may be in 
the same row or column is called an assignment problem. A set of matrix 
positions no two in the same row or column is called an independent set. 
An independent set of k matrix positions will also be called a fc-assignment. 

A random assignment problem, or RAP for short, is given by a number 
k, and an m by n matrix (min(m, n) > k) of random variables. If P is a 
random assignment problem, we denote by E{P) the expected value of the 
minimal sum of an independent set of k matrix elements. 

In this article we prove the following. 

Theorem 1.1 (Parisi's Conjecture |P98j ). Let P be the RAP where k = 
m = n and the matrix entries are independent exponential random variables 
with intensity 1. Then 

£(P) = l4 + i + - + ^. 
We also prove the following two generalizations. 

Theorem 1.2 (Conjectured by D. Coppersmith and G. B. Sorkin ICSOSj ). 

Let P be an RAP where the matrix entries are independent exponential ran- 
dom variables with intensity 1. Then 



(1) E{p) = y: 



,+,<fc {m-i){n-j) 
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Theorem 1.3 (Conjectured in |LW00j ). Let P be an RAP where some 
matrix entries are zero and all the other entries are independent exponential 
random variables with intensity 1. Then 

(2) E{P) = —J2 ' 

where di,j{P) is an integer coefficient defined in terms of the combinatorics 
of the set of zeros, see Section [KR 

The identity ^ will be referred to as the cover formula. When P has no 
zeros, 

Hence Theorem \1.'A\ =^ Theorem 11.21 To see the implication Theorem 11.21 ^ 
Theorem 11.11 (proved in |(yS98j ) , note that if we put m = n = k in then 
the terms for which gcd(A: — i,k — j) = d sum to l/d^ for d = 1, . . . ,k. 

As a consequence of Theorem 11.11 we also obtain a new and completely 
different proof of the following theorem, conjectured by M. Mezard and 
G. Parisi |M P85j . 

Theorem 1.4 (D. Aldous [IMlEIn] ). Let P be the RAP where k = m = n 
and the matrix entries are independent exponential random variables with 
intensity 1. Then 

lim E{P) = — . 

k--*oo 6 

In Section d we mention some other corollaries of Theorem 11.81 

1.1. Outhne of the proof. A key result, proved in Sectional is a formula 
for the probability that a nonzero row in an RAP is used in the optimal 
A;-assignment, see Theorem 14.11 Prom |LW00j we know that the probability 
that a nonzero element in an RAP P is used in the optimal /^-assignment 
can be written E{P) — E(P'), where P' is obtained from P by setting the 
matrix element in question to zero, see Theorem 12. 21 Therefore, the formula 
for the probability that a row (or column) is used gives linear recursions for 
the values of the RAP's. Provided m or n is sufficiently large compared to 
k, this system of linear recursions has a unique solution given by the cover 
formula in Theorem 11.31 see Section |SJ Finally in Section we prove that 
for fixed k, fixed m, and a fixed set of zeros, E(P) is given by a rational 
function in n, which must then agree with the cover formula. 

1.2. Background. Random assignment problems have attracted the atten- 
tion of researchers from physics, optimization, and probability. There are 



experimental results in |U921 lPR93j . Constructive upper and lower bounds 
on E{P) h ave bee n given in |W79l 10921 IHKMl lL93l lK87l Kj;K93| . Mezard 
and Parisi |MP85j used the non-rigorous replica method and arrived at the 
conjectured limit tt^/G. This limit was subsequently established rigorously 
by Aldous |A92l lAOlj using the weak convergence method on a weighted 
infinite tree model. In this paper we continue the exact formulas-approach 
inspired by |P98] and developed further in (ASn2. .BCeM IHKMl EHOl 

iLWon, EEgoH. 



A PROOF OF PARISrS CONJECTURE 



3 



There are also interesting results on similar problems such as finding a 
minimal spanning tree in a graph with random edge weights BFM98,, FM89I 
IF851 lEESOlj . An intriguing question is why the Riemann (^-function appears 
in the limit of both the spanning tree and bipartite matching problems. 

One tool that we believe can be useful to a wider range of problems is 
Theorem 12.21 below and its generalization Theorem 7.3 of |LWOO) . 

2. Preliminaries 

2.1. Probabilistic preliminaries. We say that a random variable X is 
exponentially distributed with intensity a if Pr{X > t) = e~"* for t > 0. 
The intensity of X is denoted I{X). We have E{X) = l/I{X). 

A random assignment problem P is called standard if the matrix entries 
are either zero or independent exponentially distributed with intensity 1. A 
standard RAP is determined by the numbers k, m and n, and the set Z of 
zero elements. 

The following is a well-known lemma. 

Lemma 2.1. Let Xi,...Xn be independent exponential variables with in- 
tensities ai,...,an respectively. Then the probability that Xi is minimal 
among Xi, . . . Xn is 

ai + • • • + a„ 

The minimum is an exponential variable Y of intensity ai + • • • + a„ which 
is independent of which variable is minimal. Under the condition that Xi 
is minimal, Xi = Y , and for j ^ i, we can write Xj = Y + X'-, where X'- 
is exponential of intensity Ui, and the variables Y and Xj for j ^ i are all 
independent. 

We say that a deterministic assignment problem is generic if no sum of 
a set of nonzero matrix elements is equal to the sum of a different set of 
nonzero elements. In the RAP's that we consider, the distributions of the 
nonzero matrix elements are continuous. Hence all RAP's considered here 
will be generic with probability 1. In the generic case, a nonzero element is 
used either in every optimal /c-assignment, or in none. Hence we can without 
ambiguity speak of the probability that a certain nonzero element is used in 
the optimal /c- assignment, without specifying whether we take this to mean 
some optimal fc-assignment, or every optimal /c-assignment. 

Notice that this is not the case for zero elements. Even in a generic 
case, there may be several different optimal fc-assignments, that differ in the 
choice of zero elements. 

The following theorem, Theorem 2.10 of |LWf)nj . is essential for the re- 
cursion equations in Sectional 

Theorem 2.2. Suppose P is a standard RAP. Suppose that the entry in 
position {i,j) is not zero. Let P' be the standard RAP where we have replaced 
the entry in position {i,j) in P with a zero. Then the probability that {i,j) 
belongs to an optimal k-assignment in P is E{P) — E{P'). 
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2.2. Covers. We will consider sets of rows and columns in the m x n-matrix. 
A set A of rows and columns is said to cover a set of zeros Z if every matrix 
position in Z is either in a row or in a column that belongs to A. A cover 
with s rows and columns will be called an s-cover. k — 1-covers will be of 
particular importance. By a partial k — 1-cover of Z, we mean a set of rows 
and columns which is a subset of a A; — 1-cover of Z. 

By a cover of an RAP we mean a cover of its set of zeros. The cover 
coefficient dij{P) is the number of partial k — 1-covers of P with i rows and 
j columns. For this to be nonzero, i and j have to be nonnegative integers 
with i+j < k. It is convenient to regard the cover coefficient as well-defined, 
but zero, for integers i,j outside this range. 

We say that a set A of rows and columns is an optimal cover of P, if A 
covers P, and A has minimal cardinality among all covers of P. The following 
lemma is well-known. For a general introduction to matching theory we refer 



Lemma 2.3 (Lattice property of optimal covers). The set of optimal covers 
of a set of matrix positions forms a lattice, where one of the lattice operations 
is taking union of row sets and intersection of column sets, and the opposite 
lattice operation is taking intersection of row sets and union of column sets. 

In particular, there is a row-maximal optimal cover containing every row 
that belongs to some optimal cover, and similarly a column-maximal optimal 
cover containing every column that belongs to some optimal cover. 

An RAP can be reformulated in a setting of bipartite graphs with random 
weights on the edges. An assignment is then a matching. In this setting a 
cover of rows and columns is a vertex cover of the subgraph of edges with 
weight zero. We will borrow terminology from matching theory and speak 
of paths, etc. referring to the corresponding graph concepts. 

The following two theorems were together with Theorem 12.21 our main 
tools for computing E{P) recursively in LWOO . 

Theorem 2.4. Let A be a deterministic assignment problem. If there is no 
set of k independent zeros, then each row and column which belongs to an 
optimal cover of the zeros must be used in every optimal k-assignment. 

Let c be an optimal cover, and let x be a positive real number smaller 
than or equal to the minimum of the elements not covered by c. Let A' be 
obtained from A by subtracting x from the elements not covered by c, and 
adding x to the doubly covered elements. Then x has been subtracted from 
the optimal k-assignment k — \c\ times, and consequently 



The previous theorem will be used for random assignment problems, by 
conditioning on where the minimal non-covered element is. 

Theorem 2.5. Suppose P is a standard RAP. Suppose also that a column 
c is used in every k — 1-cover of P. In particular, this is the case if at least 
k elements in c are zero. Let P' be the RAP with column c deleted from P 
and k and n decreased by one. Then 



Similarly, if a row belongs to every k — 1-cover of P, it can be deleted. 



to [LPSfij. 



E{A) = {k-\c\)x + E{A'). 



(3) 
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3. Combinatorics of two optimal /c-assignments 

In this section we prove some results that are needed in Section [l] Con- 
sider a deterministic assignment problem. If there are two different k- 
assignments fj, and z^, their symmetric difference fiAv := U will 

be of special importance. By a (fj,,!^)- alternating path we mean a sequence 
{h,ji), (^2, J2); • • • , {irijr) of matrix positions, where positions belonging to 
jjL alternate with positions belonging to v and ix = ix+i or jx = jx+i, for all 
X = 1, . . . ,r — 1. The parameter r will be called the length of the path. It 
is easy to see that the positions in /iAz^ will form i/)-alternating paths in 
the matrix. 

If two A:-assignments fi and differ only at zero positions we say that fi 
and v are equivalent, ^ = v. 

A deterministic assignment problem is semi-generic if there is exactly 
one nontrivial sum of nonzero elements that equals a sum of a distinct set 
of nonzero elements. 

Lemma 3.1. Let A he a deterministic assignment problem. Let fi be an 
optimal k-assignment, and let a be an element which does not belong to fi. 
Suppose that there is another optimal k-assignment v that contains a. Then 
there is an optimal k-assignment v' containing a, such that fiAi^' consists 
either of one path, or of two paths of odd length. 

Proof. Let T be the path of /^Az^ that contains a. If T is of even length, then 
we let v' = fiAT. If T has odd length, then it contains one more element of 
one of fi and than of the other. Since /i and z^ have the same size, there 
must be another path T' of /xAz^ that balances, so that T UT' has equally 
many elements from /i and z^. Then we let z/' = ^A{T U T'). Clearly v' has 
the desired properties. □ 

Lemma 3.2. Let A he a semi-generic matrix with two inequivalent optimal 
k-assignments fx and v. Then fiAv has at most two paths with some nonzero 
element, and if there are two, they both have odd length. 

Proof. By Lemma 13.11 there is an optimal fc-assignment u' that contains 
a nonzero element not in fi, such that fiAv' has either one path, or two 
paths of odd length. Since there are only two equivalence classes of optimal 
^-assignments, we must have u' = u. □ 

Lemma 3.3. Let A be a deterministic assignment problem. Let ri and r2 be 
rows, and letci and C2 be columns. Suppose that A{ri,ci) = 0, A{r2,ci) > 0, 
and A{r2,C2) > 0. Suppose further that there is an optimal k-assignment 
that uses (^2,02), but not row ri. Then there is no optimal k-assignment 
that uses (r2, ci). 

Proof. Let fj, be an optimal /c-assignment that uses (r2, C2), and suppose that 
there is an optimal /c-assignment u that uses (r2 , ci ) . By Lemma 13.11 we 
may assume that /xAz/ consists of at most two paths, and that if there are 
two paths, both contain an odd number of matrix elements. 

There must be a zero element (si,ci) G /i, otherwise (r2,C2) could be 
replaced by (ri,ci) in ^u. Let S be the path of (/xAz/) \ {(r2,ci)} that 
contains (si,ci). In other words, S is the set of matrix positions in the 
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i/)-alternating path that starts at (si,ci) and continues in the direction 
opposite to that of (r2,ci). 

There must be an element {ri,d) of v since otherwise (r2,ci) could be 
replaced by (ri,ci) in i^. Since /i does not use row ri, the element {ri,d) 
is in a path of ^Az^ which is a path with one end belonging to i'. If the 
other end belongs to fi, the path has an even number of elements, and must 
therefore be the same as the path containing (r2,ci). Otherwise the other 
end too belongs to v. Then the other path has both ends in fi. In either 
case, S must end with an element of in a row which is not used by v. 

If the sum of the matrix entries in /u H S is greater than the sum of the 
matrix entries muDS, then the cost of the /c-assignment (fiAS) U {(ri, ci)} 
is smaller than the cost of /i, contradicting the optimality of /x. Otherwise 
the sum of the matrix entries in // n 5 is smaller than or equal to the sum 
of the matrix entries in n S*. Then the cost of {i^AS) \ {(^2, ci)} is smaller 
than the cost of z^, again a contradiction. □ 

Lemma 3.4. Let A be a semi-generic deterministic assignment problem, 
and suppose that there are two non- equivalent optimal k- assignments /i and 
V . Assume further that 

(1) n has a nonzero element in the last row and 

(2) u does not use the last row. 

Then there is a unique row s such that 

(1) v has a nonzero element in row s 

(2) There is a k-assignment ^' = fj, that does not use row s. 

Proof. Existence: Choose a fc-assignment fj,' = fj, which has as many matrix 
positions as possible in common with u. Since the last row is used by /i' 
but not by v, there has to be a row s which is used by v but not by /i'. If 
u has a zero element in row s, then /i' must contain a zero element in the 
same column. By replacing this zero by the zero in row s, we would obtain 
a /c-assignment /j," = ji' which has one more element in common with z/, a 
contradiction. Therefore, the element in row s which belongs to v must be 
nonzero. 

Uniqueness: Assume on the contrary that there are two different rows s 
and t that contain nonzero elements of z^, say (s,ci) and (t, C2), and two 
equivalent fc-assignments ^ such that /i does not use row s and //' does 
not use row t. If /x would not use row t either, then filXu would have two 
distinct paths both with one end in a nonzero element of z/, in the elements 
(s,ci) and (t, C2). This contradicts Lemma|3I21 

Hence there is a zero element (5,03) of /i', and similarly a zero element 
(t,C4) of ^. Note that the symmetric difference /iA/i' contains only zeros 
and therefore it must consist of a number of (/i, /x')-alternating paths, all of 
even length. Let U be the one that contains (^,04). Unless U ends at row 
s, jiAU will be a A;-assignment equivalent to avoiding both rows s and t. 
We have already seen that this is impossible. 

We may therefore assume that /iA/i' consists of a single path from (^,04) 
to (5,03). Note that for every row r which is used in this path, there is a 
A:-assignment = /i that avoids row r. This is obtained by choosing the 
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zeros from in the part of U that goes towards (t, C4), and choosing them 
from fi' in the part that goes towards (5,03). 

Let (m, d) be the position in the last row used by fi and let L be the {fi, v)- 
alternating path containing {m,d). First note that if L contains (s,ci), it 
has to end there and thus be of even length, which means that it has passed 
through (t,C2) first. 

Case 1. L intersects U. In this case, the first element of L (starting from 
(m, d)) that belongs to U must be an element of /j.. Suppose that this element 
is in row r. Then let fi" = ^ avoid row r. It follows that /x"Az/ contains 
a path of even length starting at {m,d), and another path containing the 
nonzero element (s,ci), contradicting Lemma 13.21 

Case 2. L does not intersect U. Let S be the path in /iAi/ containing 
(s,ci), {t,Ci) and (t,C2) in this order. 

We know that U and 5 intersect in (t, C4). Of the positions in U Ci S, 
let (r, c) be the one which is closest to (s, C3) in U, and let n" = fi avoid 
row r. Then ji" /\i> will consist of a cycle containing the nonzero position 
(s, ci), and two other paths, one containing (t, C2) and one containing (m, d), 
contradicting Lemma 13.21 □ 

4. The probability that a row is used in the optimal assignment 

Crucial for our proof is the following formula for the probability that a 
certain row without zeros is included in an optimal ^-assignment. Let q{P) 
denote this probability. 

Theorem 4.1 (Row Inclusion Theorem). Let P he a standard RAP, and 
let r he a row without zeros. The probability that some element in r belongs 
to the optimal k- assignment is 



where di^ is the number of partial k — 1-covers of i rows not containing the 
row r. 

To simplify notation we will assume that the nonzero row in the Row 
Inclusion Theorem is the last row in the matrix. A consequence of this 
formula is that the probability of using a certain row without zeros in an 
optimal /c-assignment does not change if further zeros are introduced in a 
row which belongs to an optimal cover. This observation turns out to be 
sufficient for the proof of the formula. 

Lemma 4.2. Let P be a standard RAP, and let r be a row that belongs to 
an optimal cover. Let P' be an RAP obtained from P by inserting another 
zero in row r. Suppose there is a row without zeros in P. Then 



Proof. Suppose that P is a standard RAP where the first row belongs to an 
optimal cover, and that there is an element in the first row, say (1, 1), which 
is not zero. Suppose further that the last row contains no zeros. We want 
to show that if we replace the element (1,1) by zero, the probability that 
the last row is used in the optimal fc-assignment does not change. 





q{P) = q{P'). 
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Let Qphe the probability space of all assignments of values to the random 
variables in P, that is, the space of all real nonnegative m by n matrices 
that have zeros in the positions where P has zeros. If A G fip, we let q{A) 
be 1 if the last row is used by an optimal fc-assignment, and otherwise. 
We let Ax, for nonnegative real x, denote the matrix obtained from A by 
setting the entry in position (1, 1) to x. 

We construct a measure preserving involution on Op with the property 
that (except possibly on a subset of probability zero) if A G rip is a matrix 
where the last row changes between being used and not being used in the 
optimal /c-assignment when A{1, 1) is set to zero, then in ^p{A), the change 
goes the other way. In other words. 



Let A G rip. If q{A) = q{Ao), we let (p{A) = A. Otherwise notice that 
if we slide the element in position (1, 1) continuously down to zero, there 
can be at most one point at which the location of the optimal /^-assignment 
changes, and at this point, the element in position (1, 1) goes from not being 
used to being used. 

At the point Af where the change occurs, the matrix is semi-generic with 
two non-equivalent optimal /c-assignments. Let /j, be the one that contains 
an element of the last row, and let be the one that doesn't. By Lemma 
11^. 4( there is a unique row s such that i' contains a nonzero element in row 
s, and so that there is a /^-assignment jj,' = ^ that contains no element of 
row s. Notice that s cannot be the first row, since the first row is used by 
every optimal /c-assignment. 

We let y?(j4) be the matrix obtained from A by swapping the elements 
in row s with the corresponding elements (elements in the same column) in 
the last row, except in the columns where row s has zeros. 

In the analysis of this mapping, we introduce some auxiliary matrices. 
Let j4j be obtained from Aj hy setting the elements in the last row that 
are in columns where row s has zeros, to zero. By Lemma 13.31 we can set 
these elements as small as we please without changing the location of the 
optimal assignments. Hence there will be no /c- assignment of smaller cost 
than ^ and v in A'j. Then since neither nor v uses any zero element in 
row s, these can of course be increased without changing the optimality of 

and V. We let A"j be the matrix where the zero positions in row s are 
changed to the values of the corresponding elements in the last row of A. If 
A is generic, every optimal /c-assignment in A"j is equivalent either to /i or 
to u. 

Now '^{A) f is obtained from A'jr by swapping row s with the last row. This 
means that in ip{A)f, there are exactly two equivalence classes of optimal 
/^-assignments, one that includes position (1, 1) and one that doesn't. But 
since the last row has been swapped with another row, we have, for x ^ f, 



13.41 but with the roles of fi and i' interchanged. It follows that ip{f{A)) = A, 
and in particular that ip is invertible. The mapping Lp is piecewise linear. 



(5) 
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and on each piece, it is a permutation of variables. Hence it is measure 
preserving. This completes the proof. □ 



Proof of Theorem \4-l\ We now establish the Row Inclusion Theorem by an 
inductive argument. Let P be a standard RAP, and suppose the last row 
contains no zero element. Notice that by Konig's theorem (see LP86j), the 
formula holds whenever P has a /c-assignment of only zeros. Suppose that 
the formula has been established for every standard RAP with fewer nonzero 
elements than P. 

Case 1: Some optimal cover of P contains a row r with at least one 
nonzero element. Then every optimal fc-assignment must use row r. Let P' 
be as P but with another zero in row r. By Theorem 14.21 the probability 
that the last row is used is the same in P as in P' . Hence we only have to 
show that difi{P) = difi{P') for i = 0, . . . , /c — 1. We have to show that if a 
set of rows can be extended to a /c — 1-cover of P, then it is possible to use 
row r in this k — 1-cover. This will follow if we can show that row r belongs 
to the optimal cover of the remaining zeros. 

It suffices to show that if a row s 7^ r is deleted from the matrix, row 
r still belongs to the optimal cover of the remaining zeros. If the deletion 
of s does not decrease the maximal number of independent zeros, this is 
obvious. Suppose therefore that the deletion of s decreases the number of 
independent zeros. Then s belongs to an optimal cover of P. Hence the row- 
maximal optimal cover of P contains both r and s. When row s is deleted, 
the remaining rows and columns including r will constitute an optimal cover 
of the remaining zeros. 

Case 2: There is a row in P with only zeros. Let P' be obtained from 
P by deleting this row and decreasing k and m by 1. The probability that 
the last row is used by an optimal assignment is clearly the same in P as 
in P' . By induction, the probability that the last row is used in P' is (with 
k = k{P),m = m{P)) 

i=o \ I ) 

We have 

for every i. Hence 



d,fl{P)=d,fl{P')+d,.i^o{P'), 



(6) 1 sr^ di,o{P) ^ 1 ^ d,,o(PO + d,_i,o(PO 



i=0 



1_ difijP') ^ 1 difiiP') 



lm—\\ 
\ i+1 ) 



2. ^.o(P ) |^^^„,_,^ + j ^ _ ^ 2. • 

Case 3: There is a unique optimal cover c consisting of only columns. 
Then these columns will be used by the optimal A;-assignment. Therefore 
the number of elements not covered by c in the optimal A;-assignment is 
independent of the random variables in the matrix. We condition on the 
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position of the minimal element not covered by c. If in each case we subtract 
this minimum from all elements not covered by c, the same nonzero elements 
will be used in the optimal ^^-assignments. Here we are using Theorem l2.4l in 
the special case of no doubly covered elements. Let P* be the RAP obtained 
by conditioning on the minimal element in P not covered by c being in row 
t, and subtracting this minimum from all elements not covered by c. Then 
P* is a standard RAP, and the new zero occurring in row t means that row 
t belongs to an optimal cover of P*. Hence in case the new zero is in the 
last row, that row must be used in the optimal A:-assignment, while if it is 
not, we can find the probability that the last row is used by induction. We 
let P\t be the RAP obtained by deleting row t from P and decreasing k by 
1. Case 1 and 2 implies that the probability that the last row is used is the 
same in P* as in P\t. It follows that the probability that the last row is 
used in the optimal fc-assignment in P is given by 

f7) 1 1 sr-sr- d~^{P\t) 

Here we have 



Y,d,,o{p\t) = {i + m+ifl{p), 



t 

since both sides are equal to the sum, taken over all t, of the number of 
partial k — 1-covers of P that use row t (and not the last row). It follows 
that (O equals 



(s\ 1 1 ^ (i + iH+i,o(P) 

^ ' m m{m - 1) ^ p-^) 



_ I I ^ di+i,o{P) _ 1 djfljP) 

i \ t+1 ) i \ I ) 

This completes the proof of Theorem 14. II 



□ 



5. Proof of the cover formula for large m or large n 



The results of the last section will now enable us to prove that the cover 
formula ((2) holds for standard RAP's whenever either m or n is sufficiently 
large compared to k. We first prove that the cover formula is consistent 
with the Row Inclusion Theorem. 

Theorem 5.1. Let P be a standard RAP where the last row contains no 
zeros, and let Pt he obtained from P by setting the element in column t of 
the last row to zero. If the cover formula ^ holds for every Pt, then it holds 
for P. 

Proof. By Theorem 12.21 the probability that the element in column t in the 
last row belongs to an optimal fc-assignment in P is equal to E{P) — E{Pt). 
Hence 

t 



Suppose that the cover formula holds for each Pt. Then 
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lE(P) = — \^ \^ ^iJK-^tJ I _ \^ 

mn ^ ^ (^~^^ (^~^^ m ^ /"m— i^ 

t i,j \ i J \ j 



In order to prove that the cover formula holds for P, it is sufficient to 
prove that 

/q\ dij{P) ^ 1 dij{Pt) sr-^ difl{P) 

V / / > /m-l\ /n~l\ „ / J / J lm-\\(n-\\ ^ / j (m-l\ ' 

i,j \ i )y j ) ^ t i,j \ i )\ j ) i \ i I 

If we write 

and fix i and j, we see that © will follow from the identity 



Here, partial covers that contain the last row will contribute l/("7^) to 

both sides. If we let dij{Pt) denote the number of partial k — 1-covers of Pt 
with i rows and j columns that do not contain the last row, it only remains 
to show that 



dijjP) d^AP) d,,3+iiP) I 1 \-^(P^ 
or equivalently 

jd~{p) + (j + i)d~^i{p) = YjkiiPt)- 

t 

Here the first term of the left hand side counts the partial k — 1-covers 
of Pt that contain column t, while the second term counts those that don't. 
Hence © holds. □ 

Next we show that the cover formula is consistent with removing a column 
which belongs to every k — 1-cover. 

Theorem 5.2. Let P he a standard RAP, and suppose that the first column 
belongs to every k — 1-cover. Let P' he the RAP obtained from P by deleting 
the first column and decreasing k and n by 1. If the cover formula holds for 
P' , then it holds for P. 

Proof By Theorem HT^ R(P) = E{P'). We have (with n = n{P) and 
k = k{P)) 

no) 1 d,AP) ^ 1 ^ d,,,.i{P') + d,AP') 

j,jr \ I / V J / i,j \ t J \ J J 

1 ^ / di^P') , d^j{P') \ _ 1 ^ d,AP') 



= E{P') = E{P). 
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n 



Theorem 5.3. If P is a standard RAP such that max(m, n) > {k — 1)^, 
then 



Proof. We show that if m > {k — 1)^, then a standard RAP which contains 
a zero in each row must have a column which belongs to every k — 1-cover. 
The theorem then follows by induction on the number of nonzero elements 
of P. It suffices to consider the case that P has exactly one zero in each 
row. Then the total number of zeros is greater than (k — 1)^. Either there 
is a column with k or more zeros, which must then belong to every k — 1- 
cover, or there is an independent set of k zeros, and in this case there is no 
k — 1-cover. 

By symmetry, the formula also holds if n > (A; — 1)^. □ 



We know from Section [S] that for fixed k, whenever m or n is large, the 
cover formula (j^J for standard RAP's holds. In order to prove that the 
formula holds for smaller values of m and n, it is therefore sufficient to show 
that if k, m and the zero positions are fixed, and we let Pn be the standard 
RAP with n columns, then there is a rational function in the variable n that 
gives the value of E{Pn) for every n which is at least as large as k and the 
number of columns with zeros. This rational function must then be equal 
to the one given by the cover formula. We prove this by induction over a 
class of RAP's which includes not only standard RAP's. 

In an exponential RAP, all matrix elements are linear combinations with 
nonnegative rational coefficients of a set Xi , . . . , Xp of independent exponen- 
tially distributed random variables. If a variable Xi in an exponential RAP 
has intensity 1, and occurs in one and only one matrix position, and this 
matrix entry is equal to Xj, then the variable is called a standard variable, 
and the matrix position where it occurs is called a standard position. 

We introduce the concept of an RAP-sequence. The idea is to treat a 
set of similar RAP's with different number of columns in a uniform way, in 
order to prove that there is a rational expression in the number of columns 
that gives the value of each RAP in the set. 

Definition 6.1. We say that a linear function /(x) = ax-|-6 in one variable 
is k-positive, iff f{x) > whenever x > k, or equivalently, if a > and 
b > —ak. 

Obviously, the sum of two or more A;-positive functions is A;-positive. 

Definition 6.2. An RAP-sequence is a sequence Pn {n > uq) of exponential 
RAP's satisfying the following: 

(1) Each Pn is an exponential RAP in a set of variables Xn^i, . . . , Xn^pi^n)- 

(2) The numbers m and k are uniform, that is m{Pi) = m{Pj) and 
k{Pi) = k(Pj) for all i and j, while Pn has n columns. 




6. Rationality of E{P) as a function of n 
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(3) In the first uq columns, the matrix elements of the various RAP's 
differ only in that the first index of the variables is changed. In 
other words, the coefficient of Xn,i in a matrix element in the first 
no columns of P„ is equal to the coefficient of Xn^^i in the same 
position in P^q. 

(4) Beyond column no, Pn has only standard elements. 

Moreover, we say that the RAP-sequence is well-behaved if 

(1) For each i such that the variables X„^j occur in the first no columns, 
there is a fc-positive linear polynomial fi{n) such that the intensity 
of Xn,i is /i(n). 

(2) Nonzero nonstandard elements occur only in columns which belong 

to the column-maximal optimal cover of the zeros, or equivalently, 
columns that intersect every maximal set of independent zeros. 

Theorem 6.3 (Rationality Theorem). Suppose P is a well-behaved RAP- 
sequence. Then there is a rational function f{x) in one variable such that 

(1) If X is a zero of the denominator of f, then x < k. 

(2) E{Pn) = f{n) for every n > no{P). 

Definition 6.4. We say that a linear combination ui of variables Xi,. . . ,Xp 
is smaller than another if for each X^, the coefficient is smaller in ui than 
in U2- We say that ui and U2 are incomparable if neither of them is smaller 
than the other, in other words, if each of them has a higher coefficient than 
the other for some variable. 

We say that a linear combination is potentially minimal in a set of linear 
combinations, if it is smaller than or incomparable with every other. 

Let c be the row-maximal optimal cover of P. We prove the rationality 
theorem by induction on a number of parameters, in the following order: 

(1) The size of the largest independent set of zeros. An RAP-sequence 

is considered simpler if it has a larger set of independent zeros. 

(2) The number of rows in c. If the number of independent zeros are 
equal, the RAP-sequence with fewer rows belonging to the row- 
maximal optimal cover is simpler. 

(3) The set of potentially minimal nonstandard elements not covered by 
c. If 1. and 2. are equal, an RAP-sequence is considered simpler if 
fewer of the nonstandard elements not covered by c are potentially 
minimal. 

(4) If 1-3 are equal, and there are two incomparable nonstandard ele- 
ments not covered by c, then an RAP-sequence is simpler if there are 
fewer variables with different coefficients in the first two (in lexico- 
graphic order, say) incomparable potentially minimal nonstandard 
elements not covered by c. 

(5) If 1-3 are equal, and there is a minimal non-covered nonstandard 
element, then an RAP-sequence is simpler if the number of variables 
occurring in this element is smaller. 

Let P be a well-behaved RAP-sequence, and suppose that the rationality 
theorem holds for every simpler well-behaved RAP-sequence with the same 
values of m and k. We may of course assume that |c| < k. Since P is 
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well-behaved, in each row all but at most k — 1 elements are standard and 
not covered by c. 

We show that E[P) can be expressed in terms of rational functions in n, 
and values of simpler well-behaved RAP-sequences. 

Case 1: There are two or more non-covered incomparable nonstandard 
elements. Let ui and U2 be the first two (in lexicographic order). We choose 
i and j such that the coefficient of Xn^i is greater in ui, and the coefficient 
of X„j is greater in U2- Let the coefficients be oi, 02, 61 and 62 so that 
ui = aiXn,i + biXnj + . . . and U2 = a2Xn,i + b2Xnj + 

Let Q and R be the RAP-sequences obtained by conditioning on (ai — 
ci2)Xn,i being smaller or greater than (62 — 6i)X„j-, respectively. 

The intensities of (oi — a2)Xn,i and (62 — bi)Xnj are /i(n)/(ai — 02) and 
/i('^)/(^2 — ^1) respectively. The probability of (oi — a2)X„,j being smaller 
than (62 — bi)Xnj is 

fi{n)/{ai - 02) 

/i(n)/(ai-a2) + /,-(n)/(62-6i)' 

and similarly, the probability of (^2 — ^i)^n.j being smaller than (ai — a2)X„,^j 

/,(n)/(62-6i) 



is 



Therefore, 

(11) E{Pn 



fi{n)/{ai-a2) + fj{n)/{b2-bi)- 
fi{n)E{Qn)/{ai - 02) + fj{n)E{Rn)/{b2 



/i(n)/(ai -a2) + /,(n)/(62-6i) 

We show that Q and R can be regarded as well-behaved RAP-sequences. If 
we condition on (ai — a2)Xn^i being smaller than (62 — bijXnj, then we can 
write 

(ai - a2)Xn,i = Yn 

and 

(62 - bi)Xn,j =Yn + Zn, 

where Yn and Z„ are independent exponentially distributed variables. The 
intensities are given by 

ai — 02 02 — oi 

and 

^<-»)-^- 

both of which are /c-positive. 

When replacing Xn^i and Xnj by the new variables Yn and Z„, only the 
nonstandard elements are affected. Either ui is smaller than U2 in Q, or 
at least the number of variables with distinct coefficients is smaller than in 
P, since ui and U2 will get the same coefficient for Yn- For any two matrix 
elements that satisfied < U4 before we conditioned on which one of Xn,i 
and Xnj is smallest, this inequality will still hold. 

Hence Q, and similarly R, are simpler than P. By induction, it follows 
that gives a rational expression for E{P), where the denominator is 
nonzero for n > k. 
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Case 2: There are no two non-covered incomparable nonstandard ele- 
ments. Then either there is among the non-covered nonstandard elements 
a minimal one, or there are no non-covered nonstandard elements. We can 
treat these slightly different cases in the same way. 

If there is a minimal non-covered nonstandard element, let aXn^i be a 
term occurring in this matrix entry. We let Sn be a set of random variables 
consisting of aXn^i and all the non-covered standard variables (if there is no 
non-covered nonstandard element, we let S„ consist only of the non-covered 
standard elements). There are at least n — k + 1 non-covered standard ele- 
ments in each non-covered row. By grouping together the standard elements 
in each row, we can write the total intensity of Sn in a uniform way as a 
sum of /c-positive terms. Therefore the intensity of the minimum Yn of the 
terms in Sn is A;-positive as a function of n. 

We condition on the minimal element in 5„. We can then replace the 
terms in Sn by new variables Yn and Zn^i where Yn — min(5'n) and Zn^i are 
the differences between the remaining terms in Sn and the minimum. 

Since all non-covered nonstandard elements contain the variable Xn,i with 
a coefficient of at least a, we can by Theorem l2 . 4l subtr act the minimum of Sn 
from every non-covered element, and add it to the doubly covered elements. 

We then get 

(12) E{Pn) = ^-J^ + jTFT ^ mmnit)), 

where Qnit) are the new RAP's obtained by conditioning on the term t 
being smallest, and performing the change of variables and subtraction of 
the minimum, and I{t) / liYn) is the corresponding probability. In order to 
write this in a uniform way for the different values of n, we group together 
the cases where the minimum occurs in a particular row and beyond column 
no. 

By permuting columns, we may assume that in those cases, the minimum 
always occurs in column + 1. The probability for each such case is (n — 
uq) / liYn). In this way, (|12() will contain the same terms for each n > no. 

If we let S'n be the subset of Sn consisting of variables occurring in the 
first no columns, and we let be the set of standard variables in column 
no + 1, we get 



(13) E{Pn) = ^-J^ + -i- ^ I(t)E{Qn{t))+'^—^ Yl nrnQnit)). 
K n) [. n) I n) ,, 



For n = no, the second sum will be empty. However, provided we can 
show that Q{t) is a well-behaved RAP-sequence simpler than P, it will follow 
by induction that the denominator of E{Q[t)) does not vanish for n < k. 
Hence the rational expression that occurs when multiplying the probability 
(n — nQ)/I{Sn) with the expression for E{Q{t)) will vanish for n = no. 
Therefore the rational expression that results will give the correct value of 
E{P) also for n = uq. 

It remains to show that Q{t) is a well-behaved RAP-sequence, and that it 
is simpler than P. We first show that Q{t) is well-behaved. We have already 
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seen that the intensities of the variables occurring in Q{t) are fc-positive. We 
therefore turn to the distribution of nonstandard elements. Possibly, there 
are some new nonstandard elements among the doubly covered elements. If 
there is no new zero among the non-covered elements, the columns of the 
doubly covered positions of course belong to the column-maximal optimal 
cover. Suppose that a new zero occurs among the non-covered elements. If 
the new zero is in a column that belongs to the column-maximal cover of 
P, then this is still the column-maximal cover of the new set of zeros. If 
the new zero is in a column that does not belong to the column-maximal 
cover, then since it is also in a row that does not belong to the row-maximal 
cover, there must be an independent set of zeros in Q{t) which is larger than 
the largest independent set of zeros in P. Therefore, the column-maximal 
cover in P, extended with the column where the new zero has occurred, 
will be an optimal cover in Q{t). Hence in any case, the columns of the 
column-maximal cover in P belong to the column-maximal cover of Q{t). 

We now show that Q{t) is simpler than P. If Q{t) contains a new zero, 
then either it gives a larger independent set of zeros, or it has to be covered 
by a column in the row-maximal cover in Q{t). In either case, Q{t) is simpler 
than P. If no new zero occurs, this must be because the minimal term in 
Sn was a term occurring in the minimal non-covered nonstandard element. 
In this case, the number of variables in this element will decrease, again 
making Q{t) simpler than P. 

Hence ()12() gives a rational expression for E[P) whose denominator is 
non- vanishing for n > k. 

This completes the proof of Theorem 16.31 

We are finally able to give a proof of Theorem 11.31 which we restate. 
Theorem 11.31 If P is a standard RAP, then 



Proof. A standard RAP with no columns can be extended to a well-behaved 
RAP-sequence Pn by inserting more columns without zeros. Hence Theorem 
16.31 shows that there is a rational function / giving the expected value of 
Pn- Since the cover coefficient dij{Pn) can be expressed uniformly as a 
polynomial in n, the cover formula gives a rational function in n which takes 
the same values as / on the infinitely many integers n > max(no, {k — 1)^). 
Hence the cover formula must agree with /, and give the value of E{P). □ 



Besides giving a new proof of Aldous' C(2)-limit theorem. Theorem 11.31 
also makes the conjectured limits of [LWODj rigorous. It follows already 
from the conjecture of Coppersmith and Sorkin, now Theorem 11.21 that if 
a,/3 > 1, then as k ^ oo, the value of the optimal /c-assignment in a [ak] 
by [/3fc]-matrix of exp(l)-variables converges to 




7. Asymptotic results 
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where A is the triangle with vertices in (0, 0), (1, 0) and (0, 1). For instance, 
when a = 1 and (3 = 2 the hmit is equal to 

^_ (log2)2 
12 2 

It also follows that the value of a standard RAP with zeros in a region 
which is scaled up with k, m and n will converge to a similar integral. In 
particular it is shown in |LW00j that the cover formula implies that the 
limit value of a standard square RAP with zeros outside an inscribed circle 
is equal to 7r^/24. 

A result presented as a conjecture in |()92j . and proved in jAOlj . states 
that in the case k = m = n with no zeros, as n — cxd, the probability that 
the smallest element in a row belongs to the optimal assignment converges 
to 1/2. We can now give an exact formula for this probability for finite k, 

1 1 
2^ 2k' 

In the case k = m, the probability that the smallest element in a particular 
row belongs to the optimal assignment is equal to the probability that the 
smallest element in the entire matrix does. For arbitrary k, m, and n, the 
probability that the smallest element in the matrix belongs to the optimal 
A;-assignment is 

^ k{k-l) _ 
2mn 

These theorems are obtained by combining the results of [LWADj and 
Theorem 11.1^1 
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